Distributed Algorithms for Maximum Link Scheduling 
under the Physical Interference Model 
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Abstract 

A fundamental problem in wireless networks is the 
maximum link scheduling problem: given a set L of 
links, compute the largest possible subset V C L 
of links that can be scheduled simultaneously with- 
out interference. This problem is particularly chal- 
lenging in the physical interference model based on 
SINR constraints (referred to as the SINR model), 
which has gained a lot of interest in recent years. 
Constant factor approximation algorithms have been 
developed for this problem, but low complexity dis- 
tributed algorithms that give the same approximation 
guarantee in the SINR model are not known. Dis- 
tributed algorithms are especially challenging in this 
model, because of its non-locality. 

In this paper, we develop a set of fast distributed 
algorithms in the SINR model, providing constant 
approximation for the maximum link scheduling 
problem under uniform power assignment. We 
find that different aspects of available technology, 
such as full/half-duplex communication, and non- 
adaptive/adaptive power control, have a significant 
impact on the performance of the algorithm; these 
issues have not been explored in the context of dis- 
tributed algorithms in the SINR model before. Our 
algorithms' running time is 0{g{L) log c m), where 
c = 1, 2, 3 for different problem instances, and g(L) 
is the "link diversity" determined by the logarithmic 
scale of a communication link length. Since g(L) is 
small and remains in a constant range in most cases, 
our algorithms serve as the first set of "sublinear" 
time distributed solution. The algorithms are ran- 
domized and crucially use physical carrier sensing 
in distributed communication steps. 



1 Introduction 

One of the most basic problems in wireless net- 
works is the Maximum Link Scheduling problem 
(MaxLSP): given a set L of links, compute the 
largest possible subset V C L of links that can 
be scheduled simultaneously without conflicts; this 
is also referred to as the one-shot scheduling [11] 
or max independent link set problem [27]. Dif- 
ferent interference models have been developed for 
capturing conflicts in wireless networks, most of 
which are based on the idea of "conflict graphs" [24]. 
MaxLSP is challenging under most of these mod- 
els — the decision version of this problem is NP- 
Complete under many models, and constant factor 
approximation algorithms are known for many inter- 
ference models [24]. In recent years, a more real- 
istic interference model based on SINR constraints 
(henceforth referred to as the SINR model) [11] has 
gained a lot of interest. This model is much harder to 
analyze than models based on conflict graphs [3]. 

Since link scheduling is a common subroutine in 
many other problems, distributed algorithms with 
low complexity are crucial. Efficient (i.e., polylog- 
arithmic time) distributed algorithms have been de- 
veloped for many graph based interference models 
(genetically referred to as the "Radio Broadcast Net- 
work (RBN)" model) for fundamental problems, in- 
cluding independent sets, coloring and dominating 
set, e.g., [20,26]. In the RBN model, a receiver faces 
a collision if transmissions from multiple senders 
reach it simultaneously. Algorithms in this model 
rely on having spatial separation between senders to 
ensure interference free communication. 

In this paper, we focus on distributed algorithms 
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for MaxLSP in the SINR model. Even central- 
ized algorithms for this problem are much harder 
in the SINR model, than in the disk based interfer- 
ence model; recent work by [10, 14-16, 27] gives 
constant factor approximation algorithms for vari- 
ous instances of MaxLSP in the SINR model. The 
centralized algorithms of [10, 15, 27] are based on 
a greedy ordering of the links, which requires es- 
timating the "affectance," (which, informally, is a 
measure of interference), at each stage (this is dis- 
cussed formally later in Section 3) — this is one of 
the challenges in distributed solutions to MaxLSP. 
We note that efficient time distributed algorithms for 
scheduling all the links (i.e., the coloring version) is 
already known [13]. Adapting them would imme- 
diately yield a distributed O (log m) -approximation 
to MaxLSP, but it is not clear how to obtain a 
distributed 0(l)-approximation. Further, an impor- 
tant aspect of MaxLSP is that the senders and re- 
ceivers of all the links should know whether they 
have been chosen, since this is an important require- 
ment in many networking applications; this seems to 
be difficult to ensure through random access based 
approaches. 

Our contributions. In this paper, we develop a set of 
fast distributed 0(l)-approximation for MaxLSP, 
in which all nodes are constrained to use uniform 
power levels for data transmission (we refer to this 
as MaxLSP u ), improving upon the results implied 
by [2, 8]. Our algorithms and the proofs build on 
ideas from [10, 15,27] and [25], and one of the key 
technical contributions of our work is the notion of 
a "ruling", and its computation in the SINR model 
(discussed below). Our results raise two new is- 
sues in the context of distributed algorithms in the 
SINR model — adaptive power control (i.e., the fea- 
ture of using lower than the maximum power level, 
as needed), and full/half duplex communication (i.e., 
whether nodes can transmit and receive simultane- 
ously). We find these features impact the perfor- 
mance of the algorithms quite a bit. We summarize 
some of the key aspects of the results and main chal- 
lenges below. 

(1) In the case of "non-adaptive power control" 
(i.e., if all nodes are required to use fixed uniform 
power levels), we design a distributed algorithm that 
provably runs in 0(g(L) log 3 m) time and gives an 
O(l) approximation to the optimum solution for half 



duplex communication, and we improve the running 
time to 0(g(L) log 2 m) for the case of full duplex 
communication; here g(L) denotes the "link diver- 
sity", which is the logarithm of the ratio of the largest 
to the smallest link length (this is defined formally 
in Section 3). If nodes are capable of "adaptive 
power control" (i.e., they can use varying power lev- 
els for scheduling, but not data transmission), we 
improve the running time of the above algorithm to 
0(g(L) log 2 m) time for half duplex communica- 
tion, and 0(g(L)logm) time for full duplex com- 
munication. Note that in the adaptive power control 
case, the algorithm uses varying power levels during 
its run, but the links which are selected finally use 
the fixed uniform power level for data transmission. 

(2) One of our key ideas is the parallelization of 
the centralized link selection step of [10, 15, 27], 
which requires sorting all links, processing a larger 
set of links simultaneously, and efficient filtering 
based on spatial and interference constraints in par- 
allel. Moreover, it turns out that the usual notion of 
independence based on spatially separated nodes is 
inadequate because of the spatial separation of the 
sender and the receiver of a link: it is the senders 
which makes the distributed decision of transmission 
and the participation in the independent set, while 
the SINR model is receiver-oriented and it is hard 
for each sender of a candidate link to determine the 
interference caused by the chosen links at the corre- 
sponding receiver. One of the important steps of our 
algorithm involves the distributed construction of a 
"ruling" (first introduced in [7]) which relates to the 
notion of independence and aids the solution to MIS 
and coloring problems in graph topologies [4, 22]. 
The extension of the notion of ruling and its com- 
putation in the SINR model is one of the important 
technical contributions of our paper. We believe this 
basic construct would be useful in other link and 
topology control problems. 

(3) Our algorithms crucially use physical carrier 
sensing and the distributed decisions are made based 
on the Received Signal Strength Indication (RSSI) 
— this approach was used by Scheideler et al. [25] 
for distributed constant density dominating set con- 
struction. This aspect is another distinction from dis- 
tributed algorithms in the RBN model. Further, our 
algorithm uses constant size messages, and all the 
steps can be implemented within the model without 
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any additional capabilities or assumptions (e.g., such 
as those made in [2]). 

Organization. We discuss the network model and 
relevant definitions in Section 3. We present the 
high-level distributed algorithm in Algorithm 1 with 
a constant approximation ratio in Section 4. We 
introduce and analyze the distributed algorithm to 
compute a ruling in Section 5. In Section 6 we show 
the detailed implementation for each step of the high- 
level Algorithm 1; we present a second method to 
implement Algorithm 1 in Section 7, improving the 
running time by a logarithmic factor. 

2 Related Work 

There has been a lot of research on link schedul- 
ing and various related problems, because of their 
fundamental nature. Two broad versions of these 
problems are: scheduling the largest possible set of 
links from a given set (maximum independent set), 
and constructing the smallest schedule for all the 
links (minimum length schedule). These problems 
are well understood in graph based interference mod- 
els and efficient approximation algorithms are known 
for many versions; see, e.g., [24]. Distributed algo- 
rithms are also known for node and link scheduling 
(and many related problems) in the radio broadcast 
model [1, 18-20,26]. These algorithms are typically 
randomized and based on Luby's algorithm [18], and 
run in synchronous polylogarithmic time. There are 
varying assumptions on the kind of information and 
resources needed by individual nodes. For instance, 
[18] require node degrees at each step (which might 
vary, as nodes become inactive). Moscibroda et al. 
[20] develop algorithms that do not require the de- 
gree information, and run in 0(log 2 n) time. In re- 
cent work, Afek et al. [1] develop a distributed al- 
gorithm for the maximal independent set problem, 
which only requires the an estimate of the total num- 
ber of nodes, but not degrees. 

Link scheduling in the SINR model is consid- 
erably harder than in graph based models. Sev- 
eral papers developed 0(g(L))-approxfmations for 
MaxLSP u , e.g., [6, 11], which have been improved 
to constant factor approximations by [10, 15,27] for 
uniform power assignments. Some of these papers 
use "capacity" [10, 15] to refer to the maximum link 



scheduling; however, we prefer to avoid the term 
capacity in order to avoid confusion with the total 
throughput in a network, which has been tradition- 
ally referred to as the capacity (e.g., [12]). Recently, 
Halldorsson and Mitra [14] extend the 0(1) approx. 
ratio to a wide range of oblivious power assignments 
for both uni- and bi-directional links (including uni- 
form, mean and linear power assignments). This has 
been improved by Kesselheim [16], who developed 
the first 0(l)-algorifhm for MaxLSP with power 
control and an thus an 0(log ?n) -algorithm for the 
minimum length schedule problem. Most of the re- 
sults except those using uniform power assignments, 
assume unlimited power values; otherwise the results 
may degrade by a factor depending on the ratio of the 
maximum and minimum transmission power values. 

Most of the above algorithms for scheduling in 
the SINR model are centralized and it is not clear 
how to implement them in a distributed manner effi- 
ciently. The closest results to ours are by Asgeirsson 
and Mitra [2] and Dinitz [8], using game theoretic 
approaches; the former obtains a constant approx. 
ratio improving over the latter's 0{d^ ax ) approx- 
imation for MaxLSP u . Their running time can 
be much higher than ours, and they require addi- 
tional assumptions (such as acknowledgements with- 
out any cost), which might be difficult to realize in 
the SINR model. 

For the minimum length schedule problem 
(MinLSP) (where one seeks a shortest schedule to 
have all the links in L transmit successfully) un- 
der a length-monotone sub-linear power assignment, 
Fanghanel et al. [9] develop a distributed algorithm 
with an approximate ratio of 0(g(L)) times a log- 
arithmic factor. Recently, Kesselheim and Vocking 
[17] propose an O (log 2 m) -approximate algorithm 
for any fixed length-monotone and sub-linear power 
assignment. The approx. ratio of that algorithm has 
been improved to O (log to) (matching the best per- 
formance of known centralized algorithms) by the 
analysis of Halldorsson and Mitra [13], who also 
prove that if all links uses the same randomized strat- 
egy, there exists a lower-bound of O, (log to) on the 
approx. ratio. However, it is not clear how to use 
these results for MinLSP to get a constant factor ap- 
proximation for MaxLSP, in which the senders and 
receivers of all links know their status. 
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3 Preliminaries and Definitions 



Table 1: Notation. 



G 


network graph 


d{u, v) 


dist. of u and v 


V 


set of nodes 


L 


set of links 


n 


#nodes 


9{L) 


link diversity 


m 


#links 


OPT() 


optimum instance 


a 


path-loss exponent 


x(l) 


sender of link I 


P 


SINR threshold 


r(Q 


receiver of link I 


N 


background noise 


d(0 


length of link I 


A 


affectance 


SP 


sensed power 



We let V denote a set of tranceivers (henceforth, 
referred to as nodes) in the Euclidean plane. We 
assume L is a set of links with end-points in V, 
which form the set of communication requests for 
the maximum link scheduling problem at any given 
time, and \L\ = m. Links are directed, and for link 
I = (x(l),r(l)), x{l) and r(l) denote the transmit- 
ter (or sender) and receiver respectively. For a link 
set L', let X(L') denote the set of senders of links 
in V '. Let d(u, v) denote the Euclidean distance be- 
tween nodes u, v. For link I, let d(l) = d(x(l),r(l)) 
denote its link length. For links 1,1', let d(l',l) = 
d(x(l'),r(l)). Let d m i n and d max denote the small- 
est and the largest transmission link lengths respec- 
tively. Let B(v,d) denote the ball centered at node v 
with a radius of d. Each transmitter x(l) uses power 
P(l) for transmission on link I; we assume com- 
monly used path loss models [6, 11], in which the 
transmission on link I is possible only if: 



P(i) 



N(l + 



(i) 



where a > 2 is the "path-loss exponent", f3 > 1 
is the minimum SINR required for successful recep- 
tion, N is the background noise, and <j) > is a con- 
stant (note that a, f3, (p and N are all constants). 

We partition the set of transmission links into non- 
overlapping link classes. We define link diversity 
g(L) = flog 2 f^]. Partition L = {Li},i = 
1,2, . . .,g(L), where each Lj = {I | 2 % ~ x d min < 
d(l) < 2 l d m i n } is the set of links of roughly sim- 
ilar lengths. Let = 2 l d m i n , such that di is an 
upperbound of link length of Li\ and Mi, \/l £ Li, 
we define d{l) = dj. In a distributed environment, 
nodes use their shared estimates of minimum and 



maximum possible link length to replace d m i n and 
dmax, as stated in the previous section. g(L) in most 
cases < 6 log 10 and remains a constant* ; further, as 
discussed earlier, each link can compute which link 
class it belongs to. The reverse link of a link I, de- 
noted by I , is the same link with transmission direc- 
tion inverted. For a link set L' , We use L' to denote 
the set of reverse links of L' . 

Wireless Interference. We use physical interference 
model based on geometric SINR constraints (hence- 
forth referred to as the SINR model), where a subset 
V C L of links can make successful transmission 
simultaneously if and only if the following condition 
holds for each link I G L'\ 



P(l) 



2^i'eL'\{i} d<*(i',i) 



N 



(2) 



Such a set L' is said to be independent in the context. 
The Maximum Link Scheduling Problem 

(MaxLSP). Given a set of communication requests 
(links) L, the goal of the MaxLSP problem is to 
find a maximum independent subset of links that can 
be scheduled simultaneously in the SINR model. 
MaxLSP u is an instance of MaxLSP where links 
in a solution use a uniform power level for data 
transmission; note that this does not necessarily 
restrict scheduling to uniform power. In this paper, 
we use OPT(L) to denote an optimum solution 
to the MaxLSP u , and thus \OPT(L)\ is the 
cardinality of the largest such independent set. As 
discussed earlier, computing OPT(L) is NP-hard, 
and we focus on approximation algorithms. We 
say an algorithm gives a C-approximation factor if 
it constructs an independent link set L' C L with 
\L'\ > \OPT(L)\/C. 

Distributed Computing Model in the SINR-based 
Model. Traditionally, distributed algorithms for 
wireless networks have been studied in the radio 
broadcast model [5, 20, 26] and its variants. The 



* The minimum link length is constrained by the device 
dimension, empirically at least 0.1 meter; the maximum link 
length depends on the type of the network, and is usually 
bounded by 10 meters. For example, the Wi-Fi transmission 
range is below hundred meters and even long-distance Wi-Fi 
networks [23] have an experimental limit of hundred kilome- 
ters; in cellular networks the coverage is at most tens of kilome- 
ters; the transmission range in Bluetooth or 60GHz networks is 
smaller. This implies often g(L) < log 10 6 . 
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SINR based computing model is relatively recent, 
and has not been studied that extensively. There- 
fore, we summarize the main aspects and assump- 
tions underlying this model: (1) The network is syn- 
chronized and for simplicity we assume all slots have 
the same length. (2) All nodes have a common esti- 
mate of to, the number of links, within a polynomial 
factor; (3) For each link I € L, x(l) and r(l) have 
an estimate of d(l), but they do not need to know 
the coordinates or the direction in which the link is 
oriented; (4) All nodes share a common estimate of 
drain and d max , the minimum and maximum possi- 
ble link lengths; (5) We assume nodes have physical 
canier sensing capability and can detect if the sensed 
signal exceeds a threshold. As discussed in [25], this 
can be done using the RSSI measurement possible 
through the Clear Channel Assessment capability in 
the 802.11 standard. Given a threshold Thres, we 
assume that a node is able to detect if the sensed 
power strength is > Thres. 

Sensed Power-strength and Affectance. For ease 
of analysis based on links, we define affectance^ as 
that in [14, 15]: the affectance, caused by the sender 
of link I' to the receiver of link /, is A(l',l) = 

8 d a (l) 
--77T — — . Likewise, we have affectance 

i _ rf "(0 d a (l'l) 

from a set V of links, as A(L', I) = YIv&l 1 0- 
It can be verified that Inequality (2) is equivalent to 
A(L' \ {I}, I) < 1, signifying the success of data 
transmission on I. 

To simplify the analysis based on nodes, we de- 
fine sensed power-strength SP(w,v), as the sig- 
nal power that node v receives when only w is 
transmitting (which includes background noise); that 
is, SP(w,v) = P/d a (w,v) + N. Likewise, 
we have SP from a node set W: SP(W, v) = 
T,wew P/d a {w, v) + N. Let Thres(d) = P/d a + 
N be a function of distance d, such that for a node v, 
if any other node is transmitting in a range of d, its 
sensed power will exceed Thres(d). 
Node Capabilities for Distributed Scheduling 
(1) Half/full Duplex Communication: Wireless ra- 
dios are generally considered half duplex, i.e., with a 
single radio they can either transmit or receive/sense 



' Sometimes it shares the same definition with the term rel- 
ative interference, e.g., in [27]; however, "relative interference" 
may refer to other forms, e.g., in [15]. 



but not both at the same time. Full duplex radios, 
which are becoming reality, enable wireless radios to 
perform transmission and reception/sensing simul- 
taneously. (2) Non-adaptive/adaptive Power. Al- 
though links in a solution to MaxLSP u use a uni- 
form power level for data transmission, they are usu- 
ally capable of using adaptive power which vary 
across different power levels that may be used for 
scheduling. The capabilities can play a vital role in 
distributed computation. 



4 Distributed Algorithm: Overview 

In this section, we present the distributed algorithm 
for MaxLSP u . Because the algorithm is quite com- 
plicated, we briefly summarize the sequential algo- 
rithm of [10, 15,27] below, and then give a high-level 
description of the distributed algorithm and its anal- 
ysis, without the implementation details of the indi- 
vidual steps in the SINR model. Section 5 describes 
the algorithm for computing a ruling in the full and 
half duplex models. The complete distributed im- 
plementation and other details are discussed in Sec- 
tions 6 and 7, for the non-adaptive and adaptive 
power control settings, respectively. 

4.1 The Centralized Algorithm 

We discuss the centralized algorithm adapted from 
[10, 15,27] for MaxLSP u , which forms the basis for 
our distributed algorithm. The algorithm processes 
links in non-decreasing order of length. Let L be 
the initial set of links, and S the set of links already 
chosen (which is empty initially). Each iteration in- 
volves the following steps: 

(1) picking the shortest link I in L\S and removing 
I from L, 

(2) removing from L all the links in {/' € L \ S : 
A(S,l') > Co} where Co < 1 is a constant, i.e., 
all the links in L \ S that suffer from high inter- 
ference caused by all chosen links in S, and 

(3) removing from L all the links in {/' G L \ S : 
d(l', I) = c\d(l)} where c\ is a constant, i.e., all 
the nearby links of Zin L \ S. 

The results of [10, 15,27] show that: S is feasible 
(i.e., the SINR constraints are satisfied at every link), 
and | S | is within a constant factor of the optimum. 
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Consider a link / that is added to S in iteration i. The 
proof of feasibility of set S involves showing that for 
this link I, the affectance due to the links added to 
S after iteration i is at most 1 — cq, so that simul- 
taneous transmission by all the links in S does not 
cause high interference for I. The approximation fac- 
tor involves the following two ideas: (1) for any link 
I G S, there can be at most O(l) links in OPT(L) 
which are within distance cod(l), and (2) in the set of 
links removed in step 2 due to the affectance from S, 
there can be at most 0(1) links in OPT(L). We note 
that the second and third steps are reversed in [27], 
while [15] does not use the third step. However, we 
find it necessary for our distributed algorithm, which 
uses the natural approach of considering all the links 
in a given length class simultaneously (instead of se- 
quentially). Our analysis builds on these ideas, and 
property (1) holds for our case without any changes. 
However, property (2) is more challenging to ana- 
lyze, since many links are added in parallel. Another 
complication is that the distributed implementation 
has to be done from the senders' perspective, so that 
the above steps become more involved. 

4.2 Additional Definitions 
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Figure 1: Example of an (uji, LJ2) -ruling: W = {vi} U 
{ui} is the set of all dots (open and dark), while 
Ru)-l,ui(W) = {vi} which is the set of all the open dots 
denotes a (u>i, 0J2) -ruling of W. Note that all the nodes in 
W are uncovered by Rwi,u 2 (W), while all the nodes in 
-Rwi,w 2 (W) are u)\ away from each other. 



Cover and Ruling. Let W, W denote two node sets. 
We say a node u is uj-covered by W', if and only if 



3u' £ W', d(u, u') < uj; based on that, we say W is 
oj-covered by W', or equivalently W' co-covers W, if 
and only if every node in W is u; -covered by W. An 
(^1,^2) -ruling (where ui\ < 0J2) of W , introduced 
in [7], is a node set denoted by R ul>C j 2 (W), such that 

(1) VW^; 

(2) all the nodes in i2w 1)Wa (W) are at least u>\- 
separated; that is, Vu, u' <G R^x^^W)^ 
d(u, u') > oj\ \ and 

(3) W is UJ2 -covered by R u}lJU2 {W). 

Here, we have adopted a generalized definition by 
considering Euclidean distance rather than graph dis- 
tance. The concept of ruling has a vital role in our 
algorithm: it is used for choosing a set of spatially 
separated links and removing the nearby links of the 
chosen links. Figure 1 gives an example to illustrate 
these notions. 

4.3 High-level Description of the Dis- 
tributed Algorithm 

We have provided a detailed discussion of the cen- 
tralized algorithms and basic ideas in Section 4.1; 
now we discuss the distributed algorithm at a high 
level (Algorithm 1), and prove the main properties. 
We use the following constants in the algorithm 71 = 

(l^Jlf^ 1 ^) ^ + 2 and 72 as an arbitrary con- 
stant > 71, where a, /3, (f> are constants in described 
in Inequality (1) and i/j is a constant that can take any 
value from (0, 1). The algorithm sweeps through the 
link classes in g(L) phases. In the ith phase, where 
i € [l,g(L)], it selects a subset of links from Lj 
to include in S, and removes a subset of links from 
Uj>jLj to speed up later phases; the comments in 
Algorithm 1 explain each step. Step 1 in line 5-8 
of the algorithm eliminate links which do not satisfy 
affectance constraints — its implementation is for- 
mally described in Sections 6 and 7, and crucially 
relies on the measurement of the received power at 
each sender. Step 2 in line 9 of the algorithm con- 
structs the distributed ruling, which is discussed in 
Section 5 for the non-adaptive power case, and ex- 
tended to the adaptive case in Section 7. 

Lemma 4.1 (Correctness). The algorithm results in 
an independent set S. 

Lemma 4.1 shows the correctness and is proved in 
Appendix A. 1 . Algorithm 1 results in a constant ap- 
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Algorithm 1: Distributed Maximum Link 
Scheduling 
input : Set L of links 



proximation ratio shown in Theorem 4.2; its proof is 
in Appendix A.2, where this theorem is backed by 
Lemma A.2 and A.3 in Appendix A.2. The two lem- 
mas, independent of Algorithm 1, stands on its own 
to provide insights of how an optimum solution is 
shaped under the SINR model; they can be proved by 



using a combination of techniques found in [15,27]. 
We provide their proofs in Appendix A.2 to make the 
paper complete and help the reading. 

Theorem 4.2 (Approximation Ratio). V71 , 72, ip > 

0, \OPT(L)\ < C 3 ( 72 ,nS\ if 72 > 71 > 

1, where C 3 (72,V>) = ^1(72) + C 2 (Y>(1 - 
(^) 1/a n ^ere C x {x) = and 

C 2 {x) = ( ^fyja_i ) a l x + 1- are functions with 
constant output values for constant input arguments. 

5 Distributed Algorithm: (wi,^)- 
Ruling 

In this section, we present Algorithm 2, the dis- 
tributed algorithm to compute an (uj\, 002) -ruling, for 
full duplex communication under the physical inter- 
ference model; in the end of the section, we extend it 
to the half duplex setting (where a node can perform 
transmission and reception/sensing at the same time) 
with added running time. While this algorithm can 
be interesting by itself, it serves as a significant build- 
ing block for our distributed implementation. For the 
algorithm to function properly, we require the input 
parameter oj 2 > (36§5±) a ~ 2 o;i- Recall that B(v,d) 
denotes the ball centered at v with a radius of d. Let 
n be the total number of nodes. The last input pa- 
rameter b max denote the estimate of the maximum 
number of nodes in the ball B(v,uji) of any node 
v G W\ ; in the worst case, b max < n. 

In this algorithm, we call an iteration of the outer 
loop (Line 1) a phase; we call an iteration of the in- 
ner loop (Line 2) a round, consisting of the coordi- 
nation step (Lines 4 through 8) and the decision step 
(Lines 9 through 12). A node v is said to be active if 
v has not joined either R or Z; otherwise, v becomes 
inactive. 

In each round, the coordination step provides a 
probabilistic mechanism for active nodes to compete 
to get in the ruling (at Line 5). Lines 6 through 8 con- 
stitute a module to resolve the issue of sensing and 
transmitting at the same time, such that two nearby 
nodes do not both enter the ruling (i.e., Lemma 5.3). 
Next, during the decision step, a subset of active 
nodes decide to join R or Z. 

In each phase, there are C4 log n rounds, such 
that we can ensure a fraction of the node population 



output: One-shot Schedule S 

1 J<-I; 

2 foreach % = 1, 2, . . . ,g(L) do /* for each 



link class in L 
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*/ 

Ji <r- Li n J, J> <- U j> i 
<- lidi,u} 2 <- l2di, 

1st step (Lines 4-6) : 



Lj n J, 



links 



J?* 

Y>(i 

</>(! 



check 

affectance constraints, s.t. 
in Si are not subject to high 
interference by links of smaller 
lengths chosen in previous phases. 
Note: we use reverse I to check 
affectance at x(l) */ 

if i > 1 then 

{I G Ji : A(U j<i S j ,'t) < 

y 13(1+4,)) ' J' * 1 \ * ' 

( £ )lMa\ Jb j J> \ jb. 

/* 2nd step (Lines 7 & 8) : check 
spatial constraints, to obtain an 
(wi,w 2 )-ruling X(JT) of X(Jf), s.t. 
the selected links (J[) are 
spatially separated and nearby 
links (J?) of similar or larger 
lengths are excluded */ 

construct link set J\, s.t. X(Jf) is an 

(wi,w 2 )-ruliag of X{Jf); 

construct link set Jf , s.t. (1) 

J*r\J? = J?\ Jl and (2) {I G JP : 

x{l) is wi-covered by X(jf)} C J? n J\ C 

{leJf: x{l) is w 2 -covered by X(J\)}\ 

I* 3rd step (Lines 9 & 10): select 
and discard links */ 

Si+- JT,J <- J\JT- /* sele ct links 

*/ _ 

J\(J?UJ*UJ?);/ 

links */ 

— UiSi, Return S; 



* discard 
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Algorithm 2: ConstructR(wi, u> 2 , Wi,W2,b max ): 
Distributed algorithm for computing an (cj%, oj 2 )- 
ruling with full duplex radios, 
input : coi , u 2 , Wi , W 2 , b max 
output: (R, Z): an (wi, oj 2 ) -ruling of W\ 

/* Each v G W± U Wjj does the following */ 

1 for i out = 0to log b max + 1 do 

2 for i i n = 1 to C4 log n do 
if v is active then 

/* Coordination Step 

(Lines 4-8) : 1 slot */ 

U(v) <r- 0; 

if v G Wi then U(v) flips to 1 w/ 
prob. ^ ; 

if = 1 then 

w transmits and senses, J(u) <— 
the power v receives in this slot; 
if I(v) > Thres{oJi) then 
U(v) <- 0; 

/* Decision Step (Lines 9-12) : 
1 slot */ 

if C/(v) = 1 then v transmits, v 

joins R; / * inactive */ 

io else 

n v senses, I{v) <— the power v 

receives in this slot; 
if I(v) > Thres(u>i) then v 



joins Z; /* inactive */ 



13 return (R, Z); 



have either joined R or Z, and we expect the max- 
imum number of active nodes in the nearby region 
of any active node to decrease by a half (proved in 
Lemma B.3 in Appendix B.l). After each phase, the 
probability for each active node to access the channel 
and compete doubles (at Line 5). After the total of 
log b rnax + 2 phases, we have Lemmas 5.2, 5.4, 5.5 
that lead to Theorem 5.1. 

Theorem 5.1 (Correctness). Algorithm 2 terminates 
in 0(logn\ogb max ) time. By the end of the al- 
gorithm: (I) R forms an (u}\,uj 2) -ruling of W\, 
and (2) Z n Wi = W\ \ R and {v £ W 2 : 
v is ui-covered by R} C Z D W 2 Q {v G W 2 : 
v is uo 2 -covered by R}, w.h.p. 



Theorem 5.1 follows directly from the lemmas be- 
low. Lemmas 5.2, 5.3 and 5.5 prove that R is an 
{oj\,u) 2 ) -ruling of W\, w.h.p. Lemmas 5.2, 5.4, 5.5 
together shows that Z complements R in W\ and 
partially in W 2 with desired properties, w.h.p. To 
help the reading flow, we defer much of the tech- 
nical content — the proof of Lemma 5.2 (which in- 
volves Lemmas B.l, B.2 and B.3) and the proof of 
Lemma 5.5 — to Appendix B. 

Lemma 5.2 (Completion). By the end of the algo- 
rithm, all nodes in W\ have joined either R or Z, 
i.e., all nodes in W\ become inactive, w.h.p. 

Lemma 5.2 implies that Z n W\ = W\ \ R. 
We say a node v G R is "good," if and only if 
d(v,v') > lui,Vv' G R and v' ^ v. In Algorithm 2, 
When a node enters R, it makes sure that there are no 
other ones entering R within a range of u)\, and it de- 
activate all the active nodes in the same range. There- 
fore, we have the following Lemmas 5.3 and 5.4. 

Lemma 5.3 (Quality of R). All nodes in R are good, 
with probability of 1. 

Lemma 5.4 (Quality of Z: Part 1). Z contains all 
the nodes in W\ U W 2 \ R that are uncovered by R, 
with probability of 1. 

Lemma 5.5 (Quality of Z: Part 2). Further, suppose 
all nodes in R are good, then all nodes in Z are u 2 - 
coveredby R, Vw 2 > (36f5±)^u;i. 

Half Duplex Communication. Now, we assume 
that nodes are in the half duplex mode, so that they 
cannot perform transmission and reception/sensing 
at the same time. In Algorithm 2, Lines 6 through 8 
make use of the full duplex capability, such that 
Lemma 5.3 is true. To account for the case of half du- 
plex communication, if we replace the one-slot deter- 
ministic full duplex mechanism (Lines 6 through 8) 
with a randomized 0(log n)-time loop — illustrated 
by the following lines of pseudo code — we have 
Lemma 5.6 for half duplex communication as the 
counterpart of Lemma 5.3 for full duplex. The cost 
incurred includes (1) the increase in the total running 
time to obtain an (ojx, oj 2 ) -ruling by O(logra), and 
(2) a weakened statement in Lemma 5.6 compared 
to Lemma 5.3. 
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In replacement of Lines 6 through 8 in Algo- 
rithm 2 for using half duplex radios. 

l for j = 1 to C5 logn do /* resolving 
half duplex communication */ 
/* in each slot */ 
if U(v) = 1 then v transmits with prob. 
1/2; 

3 if v does not transmit then 

4 v senses, I(v) 4— the power v receives in 
this slot; 

if I(v) > Thres{uJi) then U(v) <- 0; 
/* stops */ 



Lemma 5.6 (Quality of R: Half Duplex Mode). All 
nodes in R are good, w.h.p. 

Since Lemmas 5.2, 5.4 and 5.5 remain valid, we 
obtain the following theorem for the half duplex 
case. 

Theorem 5.7 (Half Duplex). There exists a mod- 
ified version of ConstructR(a;i, 6J2, W\, W 2 , b ma x) 
for the half duplex case, such that it finishes in 
0(log 2 ?nlog b max ) time and by the end of the al- 
gorithm: (1) R forms an {uj\, 0J2) -ruling of W\, 
and (2) Z C\W X = W x \ R and {v G W 2 : 
v is oj\-covered by R} C Z n W2 Q {v G W2 ■ 
v is uj 2 -covered by R}, w.h.p. 



6 Distributed Implementation with 
Non-adaptive Uniform Transmis- 
sion Power for Scheduling 

Putting everything together, we present in this sec- 
tion the distributed implementation of Algorithm 1 
when restricted to using one uniform power level for 
scheduling. 

Theorem 6.1 (Performance). Our distributed imple- 
mentation of Algorithm 1 with non-adaptive uniform 
transmission power has the following properties: 
(1) in half duplex mode, it terminates 



0(g(L) log m) time, 
(2) in full duplex mode, 
0(g(L) log 2 m) time, and 



it terminates 



in 



111 



(3) in both modes, it produces a constant- 
approximate solution to MAXLSP U . 



For the ith phase of Algorithm 1 , we present the 
distributed implementation that works even when 
there is only one fixed power level available. We as- 
sign 72 a constant value > (36^E^) a ~ 2 li' an d let 
ui = 71 di, and L02 = 72^. The distributed imple- 
mentation goes as follows. 



Algorithm 4: CheckA(y, Y', S): Distributed al- 
gorithm for checking affectance. 

input : Link sets Y,Y',S 

output: Y a = {x(l) :l€Y, A(S, I) < 



V(i 

Y b = 



( 



j£ \l/a\a 



rn a }, 



V>(i 



{x(Z) : I G Y',A(S,l) < 

) i/a n, 



Y a = {x(l) : I G Y,A(S,l) > 



\l/a\a 



Y b = {x{l) :leY',A(S,l) > 



W ~ (^y) 1/Q ) Q } 

/* in 1 time slot: */ 

1 if I € S then x(l) transmits; 

2 else if / G Y UY' then 

x(l) senses, SP{X(S),x(l)) 
x(l) receives; 
HSP(X(S),x(l)) < 



the power 



Thres[( 



0(1- 



d(l) 



4> 



-yi a d(l)) then 



^(1+4,)' 

if I G Y then x(l) joins Y a ;ifleY' 
then x(l) joins Y b ; 

i 

if I G Y then x(l) joins Y"; if I G Y' 
then x(l) joins Y b ; 



se 



8 return Y a ,Y b ,Y a ,Y b ; 



Distributed Implementation: 1st Step: With Al- 
gorithm 4, we run CheckA( Jj, Jr*, Uj^Sj) to im- 
plement the 1st step for phase i in Algorithm 1. 
V/ G Ji U J(", on Line 3 of Algorithm 4, we get 



SP(X(S),x(l)) = Zi> 



eu 



i< iS 3 d a (l',T) 



+ N. Then, 



9 



since 



SP(X(S),x(l)) 

P0-- 

<Thresy{ 



d(l) ■ 
P/(PN)' 



a/a 



~jl/a\a ' 



d(l)) 



is equivalent to 



A(Uj<iSj, I ) 



<^(l-( 







E 



cP(Z) 

1 P/(/3Ar) I'eu^iSj a {i> , i> ) 

■jl/a^a 



the sets of links whose sender nodes_are in 
Y a ,Y b ,Y E ,Y b correspond to J?,JP,J?,JP in Al- 
gorithm 1 respectively. 

Distributed Implementation: 2nd Step: Recall 
that for a link set L' , X(L') is the set of all sender 
nodes. To implement the 2nd step for phase i in Al- 
gorithm 1. we feed b max = to to Algorithm 2 and 
run ConstructR(o;i , U2, X{J£), X(J%), m). Thus, 
we obtain an (wi, L02 ) -ruling R of X( J") and Z that 
complements R in 0(log 3 m) time for half duplex 
and 0(log 2 m) time for full duplex, Then, the sets 
of links whose sender nodes are in R, Z respectively 
correspond to X(JT),X(J?) in Algorithm 1. 
Distributed Implementation: 3rd Step: The 3rd 
step of Algorithm 1 means all the links in class Lj 
and those longer links removed in the 1st step exit 
Algorithm 1 . Because our algorithm is sender based, 
the corresponding links will quit upon the decision 
of their sender nodes in the 1st and the 2nd steps. 

7 Distributed Implementation with 
Adaptive Transmission Power for 
Scheduling 

In this section, suppose we have multiple power lev- 
els at our disposal on each node-f; we present how 
this aids the distributed implementation of Algo- 
rithm 1. 



* In this paper we only study MaxLSP u where links in a 
solution use a uniform power level for data transmission; this 
does not necessarily restrict scheduling control to using uniform 
power. The general version of the problem MaxLSP that ex- 
plores power control in both scheduling and data transmission 
in a distributed setting is a hard problem and remains open. 



Theorem 7.1 (Performance). Our distributed imple- 
mentation of Algorithm 1 with adaptive transmission 
power has the following properties: 

(1) in half duplex mode, it terminates in 
0(g(L) log 2 to) time, 

(2) in full duplex mode, it terminates in 
0(g(L) log to) time, and 

(3) in both modes, it produces a constant- 
approximate solution to MaxLSP u . 

Again, note that these adaptive power levels are 
only for scheduling in the control phase; for data 
transmission in the resulting independent set S we 
still use one uniform power level. Specifically, we re- 
quire that (1) nodes have access to a set of (g (L)) 
power levels; and (2) for each i £ [l,g(L)], there 
exists a power level Pj to use such that {^) l ^ a = 
j^di, where 73 is a constant. 

We present a second method to implement the 2nd 
step of each phase in Algorithm 1 , reducing the run- 
ning time by one logarithmic factor, by (1) perform- 
ing a preprocessing to reduce b max to some constant 
C9 in O (log to) time, (2) running Algorithm 2 with 
the constant Cg in 0(log 2 m) time with half duplex 
radios and O (log to) time with full duplex radios, 
and (3) performing a postprocessing to obtain the 
sets of links required as a result of 2nd step of Algo- 
rithm 1 in one slot. We reuse the implementation for 
the 1st and the 3rd steps from the previous section. 
Due to space limit, we put the details in Appendix C. 

8 Conclusion 

In this paper, we present the first set of fast dis- 
tributed algorithms in the SINR model with a con- 
stant factor approximation guarantee for MaxLSP u . 
We extensively study the problem by account- 
ing for the cases of half/full duplex and non- 
adaptive/adaptive power availability for scheduling. 
The non-local nature of this model and the asymme- 
try between senders and receivers makes this model 
very challenging to study. Our algorithm is random- 
ized and crucially relies on physical earner sensing 
for the distributed communication steps, without any 
additional assumptions. Our main technique of dis- 
tributed computation of a ruling is likely to be useful 
in the design of other distributed algorithms in the 
SINR model. 
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A Appendix to Section 4 

A.l Proof of Lemma 4.1 

The statement of Lemma 4.1 is equivalent to that 
V/ G S,A(S \ {I}, I) < 1. Let I be an arbi- 
trary link in S, and w.l.o.g., we assume I £ Si, 
and thus I G Li. In each phase j < i, because 
IfuJ^U JJ U JJ D J?, all the links in Uj^Lj have 
been removed from J at the end of phase j. Due to 
the 2nd step, A^Sj, t) < ^(1 - (^j) 1 /")" 
First, we show that A(Uj<jSj, /) < tp. 
For any link I' G Uj<jSj, 



/3 



cF(Z) 



1 **CT) daw, I ) 
1 P/(fiN) V ' 

< A(Uj<iSj,T) < 1. 



Hence, d(x(l'),x(l)) = d(l', I ) > 
> {^-fl-d{l), implying 

rf(*',Z) _ d(x(l'),r(l)) ^ d(x(l'),x(l))-d(l) ^ 
^7T) - d(a:(I'),x(/)) - d(x(i'),a:(I)) ~ 

1 - (^i^y) 1/a - By referring to the def- 
inition of affectance — A(L)j<i5j,Z) = 
V ^(0 anH 4fi|. .<?. V\ 



( frP(Q y 

- 1 - 7U / i -> \- \ 



P/(/3iV) 



we obtain 



Next, it suffices to show that 

A(Uj>iSj\{l},l) <i-i>. 

At phase i, oj\ = 71^; the 1st step in the phase 
ensures that when link / is added to Si, any link 
I' E Uj>jLj with d(x(l'),x(l)) < oj\ will not get 
in Si. Therefore, for the set Uj>iSj \ {/}, we have: 
(1) all the nodes in X(Uj>iSj \ {I}) have a mutual 
distance of at least oo\ ; (2) the distance from any node 
in X(Uj>iSj \ {I}) to r(l) is at least uj\ — d{l); and 
(3) u)\ — d(l) > oji/2 > 0. According to Propo- 
sition A.l, by using 71 : 



36/3 q-1 l+^ U/a 9 
1— ?/> (q-2) <f> ) 



SP(X(U^S,\{Z}),r(Q) 



p7psn ) + ^- 11 is easv to veri fy 



Proposition A.l. W G V and Vt; V, i/fij a// 
?/ze nodes in V' are at least p\ away from each other, 
(2) the distance between v and any node in V' is at 
least p 2 , (3) p 2 > pi/2 > 0, then SP(V',v) < 

36(a-l) p\ P , N 

Proof. We bound the sensed power strength by par- 
titioning the plane into concentric rings all centered 
at v, each of width p2, via a similar technique to that 
in [6,21]. Let Ring{i) denote the 2th ring (where 
i = 1,2,...), which contains every node v' that sat- 
isfies ip2 < d(v',v) < (i + l)p2', let V'(i) denote 
the subset of nodes in V' that fall in Ring(i). We 
notice the following facts: (1) For any two nodes 
v,v' G V'(i), two disk centered at v, v' respectively 
with a radius of pi/2 are non-overlapping. (2) For 
any node v G V'(i), such a disk is fully contained in 
an extended ring Ring'(i) of Ring(i), with an extra 
width of pi/2 at each side of Ring(i). The area (de- 
noted by D) of each of such disks is D = tt(pi/2) 2 . 
The area (denoted by D(i)) of Ring'(i) is 

D(t) = 7T[((i + l)p2 + Pl/2) 2 - (i P2 ~ Pl/2) 2 ] 

< 3vr(2i + \)p\. 

Using \V'(i)\ < D(i)/D < 12(2i + l)pi/pf, we 
obtain 



SI(V',w)<J2\V'(i 



i=i 



DO 



+ A 



< I2(2i + 1) p| P 

-2L, p p 2 p a 

< 36(a -1) p| P | ^ 



a 



Pi P2 



□ 



A.2 Proof of Theorem 4.2 

For a node f , we define B(v,d) as the ball centered 
at v with a radius of d. With a parameter 7 > 1, 
we then define a link set B^-(l), such that for a link 
I G Li, B^(l) contains all and only the links in {Lj : 
j > i}, with their senders in the ball B(x(l), / ydi); 
in other words, B^(l) contains the links with similar 
or longer lengths, whose senders are (7^) -covered 
by x{l). For a set V C L, B^-(L') is defined as 
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Lemma A.2 (Spatial Constraint). V7 > 1, VL' C rest of the proof we assume that > (Jfr)7ft 



L, V/ € L', \OPT{Bf{l) DL')I ^ C i(7)> 



Ci(7) 



(27+l) c 
3~ 



Proo/ Let be an arbitrary link in OPT(B^(l) n 
L'). For any link k' G OPT{B}(l)C\L'), d(k' , fc) < 
d(a;(A;'),a;(/))+d(a;(Z),x(fe))+d(/fc) < (2 7 +l)d(A;). 
Therefore, 

1 > A(OPT(5|(Z) n L') \ {k}, k) 



P- 



St/ 



d°(fc) 

fc'GOPT(B|(/)nL')\{fc} d a {k\k) 



1 



d a {k) 

pJW) 



\OPT{B^l)^L')\ [2i ^l a[k) 
1 _ d°(fc) 



>/3(2 7 + l)-"|OPT(B|(0nL / )|. 

It follows that |OPT(£|(Z)nL')| < (2 7 + l) Q /73 = 
d( 7 ). □ 

Lemma A.3 (Affectance Constraint). > a«c? 
VL',L" C L, ifV n L" = and A(L', V) > V' 
for any link I G L", fton |OPT(L")| < CfeWOI-k'L 
wAe» C 2 (^') = (tIP^j W + 1- 

Proo/ If \OPT(L")\ > 0, we can express it as 
\OPT(L")\ = b\L'\ + g, such that b and 5 are non- 
negative integers and 1 < g < \L'\. We create \L'\ 
bins, each of which has a capacity of b links; we pack 
the links in OPT(L") to the bins via a first-fit sweep 
through the links in L'\ 

(1) We order the links in L' arbitrarily; let lj denote 
the jth link in L', and let Birij denote the jth 
bin. Let set L* = OPT(L") initially; then, the 
sweep proceeds in \L'\ rounds. 

(2) In the ith round (where i = 1,2, ... , \L'\), we 
pick b links in L* whose senders are the nearest 
b nodes to the sender x{U) of the ith link in L', 
and we remove those b links from L* and put 
them into Bim. 

The completion of the above packing means that 
in OPT(L"), we have b\L'\ links "near", and g G 
[1, \L'\] links "far" from the senders of links in V . If 

2(/3b) 1 /° 



b < 



"7/3, we are done. Therefore, for the 



and we show that b < C 2 (V) - 1 = ( ff)!^ )"/^ 
in this case. 

Let Z denote a "far" link in OPT(L") that is 
out of any bins. We have for any link k G 

Birii,d(x(k),x(li)) < d(x(l), x(li)), implying that 
d(x(l),x(k)) < 2d(x(l),x(li)) due to triangle in- 
equality. Since Birii U {1} C OPT(L"), we have 



1 > A(BiniJ) = (3 



fce-Binj tT(M) 

1 P7pv) 



y l2d(a:(0,2.'(ii))+ d W^ 



>/36( 



P7pv) 
1 



2d(x(0,x(Zi))M0 + l' 



That leads to that d(Z) < ( ,^ 1 2 y Q _ 1 (i(x(Z), x(Zi)). 
Then for any link k G Birii, 

d(k,l) = d(x(k),r(l)) 
<d(x(k), x(h)) + x(Z)) + d(x(l),r(l)) 

2(/36) 1 /« 



< 



-d(x(Zi),s(0). 



(/36)V« - 1 
Since Ui-Bm, U {1} C OPT(L"), we have 

1> A(u£'jj3fni,0 = /9j] 



(3) 



L'| V 

l^k&Birii d a (k,l) 



i=l 

\Bini\d a (l) 



> 



PY. *. m (bylneq.(3)) 



i=l 



> 



b PYti'eL' d a (x(i' 



d a (l) 



, 2096)1/° x„ 



1 _ rfa (Q 

1 pJTW) 



{Pb)^ a - 1 



2(06) 



The last inequality above holds because 

1 ttWT 
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< 



and 



\OPT{L")\ < {b + l)\L'\<C 2 {xl>')\L'\. □ 



Therefore, b 



We define J a = UiJf,J a = UiJ^,J b = 

U; jf , = UiJ* J r = UiJf, J* = UiJf. J^UJb 
contains all the links removed in the 1st step in Al- 
gorithm 1 due to the affectance constraints. At each 
phase i, X(J[) is an (oj i, u 2 ) -ruling of X(J£), and 
all the nodes in X(Jf) are oj 2 -covered by X(J^); 
we choose all the links in JJ to add to S, discard all 

the links in Jf U J\ (for failing affectance check), 
and also discard all the links in Jf (because of their 
proximity to the chosen links). 
We have 

L = Ui (jf U Jfu J[U Jf) 

cu^u^uu^KJD 

= J^U J*UB|(S). 

Due to Lemma A.3, |OPT(J" U < 
C 2(^(l - {j^) 1,a T)\S\\ due to Lemma A.2, 
\OPT(Bl(S))\ < Eies\OPT(B^(l))\ < 
d (72) I S\. Therefore, 

|OPT(L)| < |OPT(J^U J*) I + \OPT(B*(S))\ 
<(Ci(72) + C 2 (i,(l - (^^)^r))\S\. 

B Appendix to Section 5 

Recall that we call an iteration of the outer loop 
(Line 1) a phase of the algorithm, and an iteration 
of the inner loop (Line 2) a round; recall that B(v, d) 
denotes the ball centered at v with a radius of d. Let 
A^iv, d) denote the set of active nodes in set W\ 
that fall in the ball B(v,d) at time point t; we will 
explicitly point t out whenever we use A^ 1 (v,d). 

B.l Proof of Lemma 5.2 

The following definitions are only involved in Lem- 
mas B.l and B. 2. Let rj be a constant > Oo 9 ^) 
In one round (which corresponds to one iteration of 
the inner loop) of Algorithm 2, let U € W\ be the set 
of nodes with U{) = 1 at line 5 in the coordination 
step. We say an active node v € W\ is "lucky" in a 



round, with to being the time that the round starts, if 
and only if 

(1) veU; 

(2) U n A^{v, r]UJ\) = {v}, i.e., v has no nearby 
active nodes in U; 

(3) SP(U \ A^ 1 (v, jjwi), v) < Thres^), i.e., to- 
tal power received from faraway active nodes is 
small. 

In a round if v gets lucky, U(v) will remain 1 till the 
end of that round, and thus will elect to be included 
in R and will cause all nodes in A^iv^uJi) to get 
into Z. 



Lemma B.l. In a round % n of phase i ou t, with t\ 
being the time that the round starts, suppose that 
for each active node u € W±, \Aj. ^(u, uii)\ < 
2log &ma*- w+l at tne beginning of the round, 
then the probability for an arbitrary active node 
in W\ to be lucky in the round is at least 

2~(log b m ax -iout+3+(2r)+l) 2 ) 

Proof. For a round at phase i out , we prove the state- 
ment in the following four steps. 

(1) First, for any active node v G W\, the prob- 
ability for v to be in U is Prob(U(v) = l) = 

2~ (log bmax— iout+2) 

(2) Second, for any active node v € W\, due 
to the packing property we upper-bound the size of 
A^^v^ui) as 



k>,wi< j( ^; 2 /2)2 k-( 



2 r)log b max —iout + l 



= (2?? + l) 2 2 

Then, because logb max + 1 > i out , 
2 iogfw-w+2 > 2 . the probability for all 
nodes other than v in A t 1 (v, rjuii) to not appear in 
U (i.e.. to remain silent) is 

Prob (A^ 1 (v, tjujx) n U \ {v} = 0) 
> (1 - Prob{U{u) = 1)) 



>(l - 2~^ ogbmax ~ iout+2 ^y A ™ 1{ - v,vu ' 1 ^ 

>(1 _ 2 -( lo S b — -w+2)^(2r,+l) 2 2 1 °g b ™— 
>(!/4)(2r?+l) 2 /2 = 2 "(2r?+l) 2 _ 
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(3) Third, for any active node v € W\, we 
lower-bound the probability that v's received power 
SP(U \ J^-{v,i}U3x),v) from outside of the ball 
B(v, rj(jj\) is "low" — i.e., below Thres(uJi) — by 

(i) partitioning the plane into concentric rings via a 
similar technique to that in [6, 21], and (ii) refer- 
ring to an (rjuii, gui) -ruling and determining the ex- 
pected number of nodes in each a ring that appear in 
U, so that we can bound the power received from all 
the nodes in the rings outside of B(v , gu\). 

We partition the plane into rings all centered at 
v, each of width guj\. Let Ring(h) denote the 
hth ring, which contains every node v' that satis- 
fies hrjui < d(v',r(l)) < (h + l)r)LUi, for each 
h = 1,2,...; let Ring a (h) denote the set of ac- 
tive nodes in Ring(h). When h = 0, Ring(0) cor- 
responds to the ball B(v,gui). Let R(h) denote 
an (rjcui, r]u>\) -ruling of Ring a (h). Then by notic- 
ing that (i) Ring a (h) C U^^^A^ 1 (V, r)U\), and 

(ii) for any two nodes v', v! 6 R(h),d(v', u') > rju)\, 
we have 



^{\U r\Ring a {h)\] 
= Yl nU(v')} 

v'£Ring a (h) 

< E nu( U ')} 



(4) 



v'eRQi) 
<\R(h)\/2. 



To bound the cardinality of P(/i), we use the fol- 
lowing facts: (i) for any two nodes v', v! € R(h), 
two disk centered at v' , u' respectively with a radius 
of r)CJi/2 are non-overlapping; and (ii) For any node 
v' € Ring a (h), such a disk is fully contained in an 
extended ring Ring'(h) of Ring{h), with an extra 
width of guji/2 at each side of Ring(h). Then, by 
referring to the ratio between the areas of Ring'(h) 
and a disk, we have |P(/i)| < 8(2h + 1); Inequal- 
ity (4) yields E { \U n Ring a (h)\ } < 4(2/1 + 1). 

Therefore, v's received power from outside of the 



ball B(v, goj\) is 

K{SP<JJ\A^{v,tw),v)} 



{E E 

/i=iu'ec(/i) 



P 



d a {v',v) 



+ N 



{\U n Ring a (h)\} 



P 



h=l 

oo 



{hqui) c 



+ N 



(2/i + l) P 



— L^i ha 



< 



h=l 

12 a - 1 P 



- + N < Thres(uJi)/2. 



g a a-2uj? 
According to Markov's Inequality, we obtain 

Prob{SP{U\AY 1 l (v,uj 2 ),v) > Thres^)) < 1/2, 
thus implying that 

Prob(SP{U\A^ 1 (v,to 2 ),v) < Thres(u x )) > 1/2. 

(4) Finally, combining the above three, 
the probability that v is lucky is at least 

2-(logf>maa;-Wt+3 + (277+l) 2 ) [— ] 

Lemma B.2. In a round ii n of phase i ou t, with t\ 
being the time that the round starts, suppose that 
for each active node u E W\, \A t x (u, UJ\)\ < 
2^sb max -t out +i . tnm> j or an arbitrary active node 

v £ W\ with 1^(1;, wi)| > 2 logbma *- i ™ t , the 
probability that v becomes inactive by the end of the 
round is at least a constant Cq, where < Cq = 

2 -(3+(2r,+ l) 2 ) < L 

Proof. In a round of phase i out , a sufficient condition 
for v to be inactive by the end of the round is that 
either v or any node in A^-{v,ljx) enters R, such 
that v either enters R or Z and exits the algorithm. 
Further, that either v or any node in A^ 1 (v,lo\) gets 
lucky satisfies this condition. Therefore, the proba- 
bility for v to become inactive in the round is at least 

Prob(v' is lucky) 
>|^4j yi (u,o;i)|2~( log b ^ax-iout+3+{2v+i) 2 ) 



>2 



-(3+(2^+l) 2 



□ 
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Lemma B.3. Let Ei . denote the event that at the 

l >OUZ 

end of the phase i out , (^^(n, < b max /2 lout 
for every active node u G W\, where t\ is the time 
that the last round in i out ends. Vi ou t, Prob(Ei out ) > 
1 — iout/ n Cl , for some positive constant C-j. 

proof by induction. At the end of the phase i out = 0, 
this is trivial. Suppose that for i out = % > 1 the state- 
ment is true; we show that it still holds for i out = 
i + 1. At the beginning of phase i + 1, if we already 
have event -Ej+i, we are done; otherwise, let v de- 
note an active node such that ^^(t^cji)] > ^ff, 
where t2 is the time that the first round in i + 1 starts. 
We show Prob(E i+ i) > 1 — ^fh as below: 



n u 7 



(1) We choose a constant C4 > 



-CV+L_ fo r the 



inner loop at Line 2 of Algorithm 2. 

(2) Under the induction assumption for phase i, 
we have: for any round ij n during phase i + 1, 
with £3 being the time that this round starts, if 
\A^ X (v, J < Tp^ff, we are done; otherwise, as 
long as > l^ 1 ^, wi)| > 7p3?f, the prob- 
ability for v to turn inactive during the round is 
at least Cq according to Lemma B.2. Then, the 
probability for v to become inactive by the end of 
phase i + 1 (consisting of C4 log n rounds) is at least 
1 - (1 - C 6 ) C4l °s™ > 1 - l/n c ? +1 . 

(3) By considering both the conditional probabil- 
ity and the fact that there are at most n such nodes as 
v, 



Prob(E i+1 ] 
>Prob{E i+1 
n 

" (1 " 



Ei)Prob(Ei) 

—)>l- i+1 

n C'7 



)(1 



c 7 



□ 



At the end of phase i out = log b max +l, with t& be- 
ing the time that the last round in phase i out ends, we 
have that with probability of at least 1 - logb £j 7 x+1 , 
for every active node u G W\, [^^(n, U3\)\ < 
1/2 < 1. That means that all the nodes in W\ 
have joined either R or Z, with probability of at least 
1 



n u 7' 



— , concluding the proof of Lemma 5.2. 



B.2 Proof of Lemma 5.5 

Proof by contradiction. Suppose there exists a node 
v G Z that is not L02 -covered by R. Therefore, all 
the nodes in R are outside of the ball B(v,u)2). We 



calculate SP(R, v) and derive the conflict. We no- 
tice the following three facts: (1) all the nodes in R 
have a mutual distance of at least oj\; (2) the dis- 
tance from any sender node in R to v is > UJ2, and 
(3) 0J2 > oji/2 > 0. Due to Proposition A.l, 



SP(R,v) <36 



a 



1 



P 



.a-2 



+ N 



<36- 



n 
a 



lwf 



a 



Thres{oj\) < Thres(oJi) 



According to the condition for v to enter Z, we 
must have SP(R,v) > Thres{uj\), where lies the 
contradiction. 

B.3 Proof of Lemma 5.6 

Suppose there is "bad" node v G R such that there 
exists a node v' G R and d(v, v') < uj\. We call such 
a v' a "bad" partner of v. The only possible situation 
for v to have a "bad" partner v' is that v and u' enter 
R in the same round; otherwise, one of them should 
have been "pushed" into Z during the decision step 
of a round when the other enters R. 

Let u denote an arbitrary node in R. The necessary 
and sufficient condition for u to be bad (or to have a 
bad partner) is that, in the coordination step of the 
round when u enters R, there exists at least one ac- 
tive node v! G W\ such that, (1) d(u, u') < c*?i, (2) u 
and v! are in U, and (3) u and v! made the same ran- 
dom binary decisions all through the C5 log m slots 
of the coordination step. 

W.l.o.g., we assume u enters R at round Zj„ in 
phase i out . Let t\ denote the time that round ij n 
in phase i ou t begins. The probability that u is bad 
equals the probability for at least one other active 
node u' G u>\) {i.e., d(v,v') < uj\) to en- 

ter R at round %i n in phase i out . 



Prob(u is bad) 
< ^ Prob{U{u) 

u'&A^ 1 (u,uii) 



1 



2C 5 logm 



< 



< 



1 K^i)j 

fjC*5 ^.log bmax— iout + 2 

1 max act ive ioGVKi (iw,a;i)| 



n 



Cr, 



4l0g 6; 



max %out~ 
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Further, due to Lemma B.3, 

Prob( max IaP^tu, Wi)| < 2 l ° sbma *- i ™ t+1 

active w£Wi 1 

i m ,f — 1 1 
>1 — > 1 . 

Therefore, 
Prob(u is good) > (1 - 4r)(! " 



n min{C5,C7-l} ' 

Finally, since R contains at most n nodes, the 
probability that there are no bad nodes in R is at least 
1 - i 

X n mm{C 5 ,C 7 -l}-l " 

C Appendix to Section 7 

Recall that we assume that the adaptive transmission 
power for control has the following properties: 

(1) nodes have access to a set of @(g(L)) power lev- 
els; and 

(2) for each i G [1, g(L)], there exists a power level 
Pi to use such that (^) 1//ct = Jsdi, where 73 is 
a constant. 

We introduce a new constant 74 > (36^5^)"~ 2 7i> 
and assign 72 a constant value > 73 + 74. For the ith 
phase of Algorithm 1, let oj\ = 71 di, 0J2 = 72 du 
w 3 = 73 dj and = 74 c?, . We implement the 2nd 
step of each phase in Algorithm 1 with the following 
three substeps 

C.l Preprocessing: Constant Density Dom- 
inating Set 

Scheideler, Richa, and Santi in [25] propose a dis- 
tributed protocol to construct a constant density dom- 
inating set of nodes under uniform power assignment 
within O(logm) slots. They define Dom(W, Pt) as 
a dominating set of a node set W with transmission 
power of Pt on each node, such that W is dt-covered 
by Dom(W, Pt), where dt is the transmission range 
under P t . Then, by "constant density", they mean 
that Dom(W, Pt) is a (9(l)-approximation of the 
minimum dominating set of W, such that within the 
transmission range dt of each node in W there are 
at most a constant number Cg of nodes chosen by 
Dom(W,P t ). 



At phase i of Algorithm 1, after the 1st step of 
checking affectance, we execute the protocol on the 
node set X{Jf) with power Pj which corresponds 
to a transmission range of W3, and thus obtain a con- 
stant density dominating set Dom(X(J^), Pi) out of 
X(Jf). Dom(X(J?),Pi) has the following proper- 
ties: 

(1) Dom(X(J?),P t )eX(J?); 

(2) dominating set: all the node in X{Jf) un- 
covered by Dom{X{Jf),Pi); and 

(3) constant density: \/v G X{Jf), 1 < 
\B(v,U3) n Dom{X(J^),Pi)\ < Cg, where Cg 
is a constant. 

C.2 Construction of Ruling X(J[) 

ConstructR(wi, w 4 , Dom(X{Jf), Pi), X( j\ ), Cg) 
produces R as an {uj\, CJ4) -ruling of 
Dom(X(Jf), Pi), and Z such that 

(1) Z n X{Jf) = Dom{X{Jf), Pi) \ R; 

(2) Z n X(Jl) D {v G X(J%) : 
v is uj\ -covered by R}; and 

(3) Z n X(jf) {v G X(Jl) : 
v is u>4-covered by R}, i.e., Z is u>4-covered by 
R. 

We argue that R is an (wi, LO2) -ruling of 
X(Jf) due to the following two properties: 
(1) R C X(Jf) and any two nodes in R are 
u>i -separated, and (2) X{Jf) is c^-covered by R. 
Property (2) can be deduced from the facts below: 
(1) X(Jf) is uncovered by Dom{X(Jf), Pi) due 
to the preprocessing step, (2) Dom{X(Jf),Pi) 
is W4 -covered by R as a result of 
ConstructR(wi, w 4 , Dom{X(Jf), Pi), X{ j\ ), Cg), 
and (3) W4 + W3 < L02 by our construction. There- 
fore, R corresponds to X(Jf) in the 2nd step of 
Algorithm 1. 

C.3 Postprocessing: Accounting for X(J-) 

Construct Z' = Z U (x{Jf) \ R\ ; the following is 
true: 

(1) Z'C\X(J?) = X{J?)\R; 

(2) Z' n X(JP) 2 {v G X(jf) : 
f is -covered by R}; and 

(3) Z' n C {u G X(J^) : 
v is a>2-covered by R}, i.e., Z' is ^-covered by 
R. 
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Therefore, Z' corresponds to X(J?) in the 2nd step 
of Algorithm 1 . 
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